QCD Splitting/ Joining Functions at Finite Temperature in the 

Deep LPM Regime 

Peter Arnold and Qaglar Dogan 

Department of Physics, University of Virginia, 
Box 400714, Charlottesville, Virginia 22904, USA 
(Dated: September 2, 2008) 

Abstract 

There exist full leading-order-in-Os numerical calculations of the rates for massless quarks and 
gluons to split and join in the background of a quark-gluon plasma through hard, nearly collinear 
bremsstrahlung and inverse bremsstrahlung. In the limit of partons with very high energy E, where 
the physics is dominated by the Landau- Pomeranchuk-Migdal (LPM) effect, there are also analytic 
leading-log calculations of these rates, where the logarithm is \n.{E/T). We extend those analytic 
calculations to next-to-leading-log order. We find agreement with the full result to within roughly 
20% for > lOT, where E^ is the energy of the least energetic parton in the splitting/joining 
process. We also discuss how to account for the running of the coupling constant in the case 
that E/T is very large. Our results are also applicable to isotropic non-equilibrium plasmas if the 
plasma does not change significantly over the formation time associated with particle splitting. 



I. INTRODUCTION AND RESULTS 



When very high energy particles travel through a quark-gluon or electromagnetic plasma, 
the dominant energy loss mechanism is through hard bremsstrahlung or pair creation, similar 
to the cascading of high energy cosmic rays in the atmosphere or of a high energy particle 
in a calorimeter. It was long a problem of interest to calculate the rate for such splitting 
processes in the formal limit of very high temperature quark-gluon plasmas, where the 
running strong coupling as{T) can be treated as small [1-6]. The problem is complicated by 
the Landau- Pomeranchuk-Migdal (LPM) effect [7] (reviewed below): for very high energy 
particles, the quantum mechanical duration of the splitting process exceeds the mean free 
time between collisions, and so successive collisions cannot be treated independently. 

For the case of approximately on-shell massless particles traveling through an infinite 
medium, a complete leading-order analysis of such processes was carried out by Jeon and 
Moore [6] using the formalism of Arnold, Moore, and Yaffe (AMY) [8-10]. This analysis 
requires substantial numerical work to solve integral equations describing the LPM effect. 
Where possible, it's always nice to have analytic results in place of numerical results, because 
they are simpler to calculate and because they can facilitate comparison between different 
approaches. One case where analytic results can be found, explored in earlier literature, is 
the limit where the particle momentum p is so large compared to the plasma temperature T 
that the inverse logarithm 1/ ln(p/T) can be treated as a small number. We will refer to this 
as the deep LPM regime. In the limit of small ag, earlier authors have given analytic results 
for splitting processes to leading order in powers of this inverse logarithm.^ (In contrast, 
the work of Jeon and Moore made no assumption about the size of the logarithm, and holds 
both in and out of the deep LPM limit, to leading order in coupling Og.) The difficulty 
with leading-log results is that for practical purposes there is a huge difference, for example, 
between ln(p/T) and ln(p/47r^r) if p is of order 10-100 T. But a leading-log analysis will not 
distinguish between these two situations since ln(p/47r^T) = ln(p/T') -|- 0(1). Consider, for 
example, 5-20 GeV jets at early times in a RHIC collision, with temperatures of order 300 
MeV. As a general rule, one almost always needs to push expansions in inverse logarithms 
to next-to-leading-log order to get useful results. 

Our goal is to find analytic results for splitting processes to next-to-leading logarithmic 
(NLL) order in the deep LPM regime. In this paper, we make a first step towards that 
goal by computing the next-to-leading logarithmic result in the formal parametric limit that 
T -C p <^ T/[Q;g ln(Q;~^)] for weak coupling: that is, deep in the LPM regime, but not too 
deep. We will leave the more involved calculation of NLL results for higher energy particles 
P ~ ^/[<^s ^^("^s^^)] for future work. However, later in this paper (section VIA), we will see 
that NLL results for higher-energy particles will differ by only about 15% from the formulas 
derived here. In this paper we will also discuss leading-log results in the case of extremely 
high energy particles, where P is so large that there is significant difference between the 
running couplings as{P) and as{T). 

As we will discuss later, our results apply not just to equilibrium plasmas but more 
generally to plasmas with an isotropic distribution of particle momenta (in which case the 
"T" in our discussion refers to the typical energy scale of plasma particles) . We express our 
answers in terms of the splitting functions ^a<-*bc{P'-, xP, (1 — x)P) of Refs. [9, 10],^ which 



^ See Eqs. (19-20) of Ref. [11], which is based on the earUer work of Refs. [1, 4, 5]. 

^ We use the term "sphtting functions" to describe the functions 7 in the sphtting rate (1.1). Though 
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FIG. 1: Exact (solid line) and next-to-leading logarithmic order (dashed line) results for "fg^gg at 
leading order in coupling as for QCD with three massless flavors (A'^^ = 3). Jg^gg is plotted in 
units of a2pi/22-3/2 

function of P/T for various values of x, where P is the momentum of 
the initial high-energy particle and xP and (1 — x)P are the momenta after splitting. From top to 
bottom, the curves corresponds to x = 0.1, 0.2, 0.3, 0.4, 0.5. 



are defined so that the rate per particle of type i and momentum P for that particle to split 
into particles of type j and k and momentum fractions x and 1 — a: is^ 

= ^ 7a^6c(P; xP, (1 - X)P) [1 ± MxP)] [1 ± mi - X)P)]. (1.1) 

Here f{p) is the phase space distribution of plasma particles of a given type, and the factors 
1 ± f{p) are Bose enhancement or Fermi blocking factors. In equilibrium, f{p) = l/(e^^ =F 1) 
is the Bose or Fermi distribution associated with (massless) particles of a given type. Ua is 
the number of spin times color states for particle type a {i.e. 6 for a quark or anti-quark 
and 16 for a gluon). The prefactors in (1.1) are just convenient normalization conventions 
in the definition of 'ja^bc- 

We will give the formula for our next-to-leading-log (NLL) result in Sec. Ill, after a 
qualitative review of the form of the leading- log result. For now, we offer in Fig. 1 a 
numerical comparison of the NLL computation to the full leading-order-in-ag formula of 
Refs. [9, 10] for the case of gluon splitting g gg. Formally, both of these curves assume 
weak coupling and that P is parametrically small compared to T/[a^\n{a~^)]. Obviously, 
for any phenomenologically interesting values of as, the last assumption is rather unlikely to 
be valid for the momenta P ~ 10^ T shown at the far right of the plot. We extend the plot 



there is some connection, they should not be confused with the DGLAP splitting functions familiar from 
zero-temperature QCD. Also, in the notation of Refs. [9, 10], this would be 7^j,(-P; xP, (1 — x)P). 
^ Readers of Ref. [6] should note that we use the symbol F to denote rate per particle, whereas Ref. [6] 
denotes this by dT/dt. 




xP/T 

FIG. 2: As Fig. 1 but the horizontal axis shows p<^/T, which is xP/T for the x values used to 
label the curves, and the vertical axis has been rescaled by a factor of x'^/^. From top to bottom, 
the curves corresponds to x = 0.5, 0.4, 0.3, 0.1, 0.2. 

this far simply to show that the NLL curve is successfully approaching the more complete 
numerical calculation based on the same assumption, as it must. 

To summarize how well the NLL expansion works, it is useful to rescale the plot as in 
Fig. 2. Here, the horizontal axis shows the smallest of the two final momenta, 

j9< = min(xP, (1 - x)P), (1.2) 

instead of the initial particle momentum P, and the vertical axis has been scaled with a; in a 
way that shows the limiting small x behavior by collapsing small x curves atop each other. 
For p< ~ lOOT, the NLL result differs from the full leading-order formula by roughly 5%. 
For p< ~ lOT, the difference is roughly 20%. Pushing the expansion in inverse logarithms 
down to p</T ~ 1, where one would not expect it to be useful, we see that the NLL result 
gives the right order of magnitude but is off by a factor of roughly 2. We conclude that the 
NLL approximation for gluon splitting in an infinite medium is reasonable for p< > lOT. 

In the next section, we will set up our discussion by giving a brief qualitative review of 
the basic parametric scales associated with the LPM effect in a quark-gluon plasma. We 
also discuss the reason for the parametric assumption P -C T/[ag ln(Q;~^)] in our analysis: 
In this limit, it turns out that the momentum transfer due to scattering during a splitting 
process is small compared to the 0{T) momenta of typical plasma particles, which greatly 
simplifies the analysis. In section III, we present the formulas for our NLL result, followed 
by their derivation in section IV. In section V, we compare numerical results for processes 
besides the g ^ gg splitting presented above. 

We should clarify that throughout most of this paper, we work to leading order in powers 
of the coupling a^- In particular, the phrase "leading logarithm" will mean the leading 
logarithmic contribution at leading order in Og. It does not mean a sum of leading logarithms 
at all orders in a^. The phrase next-to- leading logarithm will be used similarly, and does 
not include any resummation of effects higher order in a^. 
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FIG. 3: A diagram contributing to bremsstrahlung from a single collision. 



In particular, our NLL results will formally assume that ln(P/T) is parametrically large 
but that asln{P/T) is parametrically small. In section VI, we discuss what to do when 
tts In (P/T) is not small. This is equivalent to a discussion of what renormalization scale 
should be used when evaluating as- In this section, we necessarily abandon the P <C 
T/[ag ln(a;~^)] restriction, considering much higher momentum scales, but we will only work 
to leading log order. 



II. QUALITATIVE REVIEW OF LEADING-LOG RESULT 



A. Basic Parametric Estimates 



Before presenting our next-to-leading-log results, it is useful to first qualitatively re- 
view the form of leading-log results. For definiteness, we'll consider the case of gluon 
bremsstrahlung. Splitting can occur when nearly collinear bremsstrahlung accompanies 
a small-angle scattering of the high energy particle off of a plasma particle, such as depicted 
by the diagram of Fig. 3 for a momentum P particle splitting into two particles of momen- 
tum k ^ xP and p ~ (1 — x)P. The intermediate solid line in this diagram is off-shell in 
energy by an amount of order 

6E = Ep + Ek ~ Ep^k 

~ 2p ^ 2k 2{p + k) 

^p\ + w? kl + ml ^ + fc_L P + 

2xP ^2(l-x)P 2P ■ 

For simplicity, in this review section we'll focus on the case x ~ 0.5 (that is, neither x nor 
1—x very small). In a typical bremsstrahlung process, the relative angle between p and k is 
the same order of magnitude as the angle of deflection in the initiating small-angle scattering 
process.^ For x ~ 0.5, then 

Ql + (masses)^ 

p ' (2.2) 

where Q± is the size of the momentum transfer (transverse to the high energy particle's 
direction of motion) in the underlying scattering process. Thermal particle masses are of 



For larger angles between p and k, there is a cancellation between the amplitudes for initial and final state 
radiation. In vacuum, bremsstrahlung can be logarithmically dominated by smaller angles between p and 
k, giving rise to a collinear logarithmic enhancement of the bremsstrahlung rate. This is nor relevant to 
the deep LPM limit, however, for reasons that will be reviewed in footnote 7. 
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order gT, and the relevant range of Q± is bounded below by the Debye screening mass of 
order gT. So we can simplify the parametric estimate (2.2) to 

5Er.^. (2.3) 

The quantum duration 1 / 6E of the Bremsstrahlung process is called the "formation" time 
of the Bremsstrahlung gluon: 

In the limit of large P, this time will become larger than the mean free time between 
collisions, and successive collisions can no longer be treated as quantum mechanically in- 
dependent for the calculation of bremsstrahlung. As far as bremsstrahlung is concerned, 
the multiple collisions during the formation time have roughly the same effect as a single 
collision with the same total momentum transfer Q± as the multiple collisions. So, on the 
one hand, the formation time is given by (2.4). On the other hand, it is also the time tcoU 
for multiple collisions to generate a total momentum transfer of size Q±, which is given by^ 

,-;,.„..(gj,„(|i).g! ,.(;?£). p.5) 

Here, n ~ T"^ is the density of particles to scatter from, and (Ji{Q±) ~ g'^/Q\ is the 
cross-section for a single Coulomb-like scattering with momentum transfer of order Qj^. 
The logarithmic enhancement to the rate is known as a Coulomb logarithm, which 
accounts for the fact that a given total momentum transfer can occur not only through 
a single collision, but also through multiple collisions which each individually transfer less 
momentum but occur more frequently.^ Self-consistently equating (2.4) and (2.5), 



This determines the total (5± and thence tform and tcou- Finally, we want the rate for 
bremsstrahlung. Every time increment tcou, we produce effectively one net collision from the 
point of view of a potential bremsstrahlung gluon. The cost of adding a gluon emission to 
a scattering process is paramctrically a factor^ of g^. So the rate for bremsstrahlung is of 
order g^/tcoii, except that we need to include final state factors 1 ± /: 

9^ 

Tbrem ~ X (final statc factors). (2.7) 

^coU 



^ For a review of these scales in extreniely simple language, see Sec. 4.5.1 of Ref. [12] in the context of 
P ~ T. For a more serious discussion in the deep LPM regime, see for example Sec. 4.1 of Ref. [2], where 
the translation to the present discussion is that their mean free path A for 2— >2 scattering is of order 
l/(p^T) and their inverse screening length fi is of order gT. 

^ For a textbook discussion of Coulomb logarithms, see Ref. [13]. 
For bremsstrahlung from a single collision in vacuum, there is also generally both an infrared logarith- 
mic enhancement, associated with integrating over small x, and a coUinear logarithmic enhancement, 
previously mentioned in footnote 4. The infrared behavior need not be considered here because we are 
restricting attention to a; ^ 0.5 in this qualitative discussion. (Note also that the discussion of the in- 
frared in medium would be different from the vacuum case because of the final state factor 1 + f for the 
bremsstrahlung particle.) The collinear logarithm in the vacuum case depends on the initial and final 
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Our qualitative discussion in this section has been for x ~ 0.5, and the Fbrcm above is 
imphcitly integrated over x's of this order. Comparing to (1.1), the sphtting functions that 
we defined earher are then of order 

7 ~ 7^ ~ g^Ql- (2.8) 

IcoU 

We can now determine Q± self-consistently from (2.6) and thence the sphtting function via 
(2.8). 



B. Preview of the inverse log expansion 



To make the leading-log expansion more explicit, it's helpful to note that the Debye mass 
mo is of order gT and rewrite the parametric equation (2.6) for Q± in the form 

gi~(0'^'lnV^(gi), (2.9) 

where 

Q±^^, (2.10) 

and then expand the solution for Q± by iteration. First ignore the logarithm and define Q±o 
by the corresponding solution 

Qlo - [j) ■ (2-11) 
Plugging this into the right-hand side of (2.9) gives an improved solution 

1/2 



Ql 

So 



(0 ln^/'(Qio) (2.12) 



and the splitting function (2.8) is 

7 - g'mlQl - (0'''lnV^ (^0 . (2.14) 



The inverse logarithm expansion is predicated on the assumption that the logarithm in 
(2.12) and (2.13) is large. As an example, supposed we plug the approximation (2.12) back 



state particles traveling in straight lines for sufficiently long distances before and after the collision. It 
does not arise in the LPM regime because t^ow ~ iform- In any case, we will shortly turn to specific 
results rather than rough parametric estimates, and the only logarithm that one finds in the answer is the 
Coulomb logarithm in (2.5). 
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into the right-hand side of (2.9) to get a new approximation for Q±: 




(2.15) 



In the hmit that ln(P/T) :$> 1, then ln(P/T) is also large compared to lnln((5io) ~ 
lnln(P/T). So we may formally choose to think of the double-log correction as part of 
the next-to-leading order terms in the expansion in inverse logs. 



C. The parametric assumption P <^ T/[a^ ^^{'^s 

As we'll discussion in section IV A, the NLL calculation will be much simpler if the 
momentum transfers in individual scattering events can be treated as a small perturbation 
to typical particles in the plasmas. The relevant transverse momentum transfers q± of 
individual collisions range from order rriD to order Qj_, which gives rise to the Coulomb 
logarithm in (2.5). So the simplifying assumption that g_L -C T is equivalent to Q± <C T in 
our analysis. Prom (2.9), this is 

which in turn is 

In weak coupling, the simplifying assumption that momentum transfers are small compared 
to T is therefore parametrically P <ti T/[ag ln(Q;^^)]. 

III. OUR NLL RESULT 

Our next-to-leading-logarithm result can be summarized in the following form: 

7,^,,(P; xP, (1 - x)P) = ml /ii(l, X, 1-x; A, A, A) ^ "^^^ , (3.1a) 

7,^,,(P; xP, (1 - x)P) = ml /ii(l, x, 1-x; F, A, F) ^^.^^Z^]' ' (^.Ib) 

7,«,,(P; xP, (1 - x)P) = ml Ulil, X, 1-x; A, F, F) (3.1c) 

(per quark flavor). 
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where iJ>]_{xi, X2, x^; si, S2, S3) solves the equation 



.2 9T 
mo 



"2 


P" 


1/2 


- X1X2X3 


T 




TT 





|(^S2 + ~ dsi)Xi In 



+ |(^S3 + C'si ~ Cs2)xl In ( J + ki^si + C'sj — Cs3)a;3ln ^^'^ 



1/2 



(3.2) 



^ is the constant 

i = exp (2 - 7e + f ) ~ 9.09916 , (3.3) 

7e is the Euler-Mascheroni constant, and the other constants will be defined in a moment. 
Schematically, these equations arc of the form (2.9) and (2.14) with fi± playing the role 
of Q±. The appearance of three different logarithms can be roughly understood as arising 
because there are three different particle momenta in the problem: P, xP, and (1 — x)P. 
(In contrast, the earlier qualitative discussion was only about orders of magnitude and took 
X ~ 0.5.) The various explicit functions of x at the ends of Eqs. (3.1) are simply the usual 
vacuum DGLAP splitting functions divided by x{l — x). The role of the medium, both in 
terms of providing the momentum transfer and the LPM effect, is contained in 7n^jj,\. 

The reason we have switched notation from Qj_ to fi± is because there are several different 
transverse momenta in the problem, associated with the three different particles, and each 
of them has a distribution of values rather than a single well-defined value. We did not want 
to give the impression that fi± exactly corresponded to a particular transverse momentum in 
the problem. For x ~ 0.5, /1± is the same order of magnitude as the Q± scale we identified in 
earlier discussion, but one could just as well redefine the normalization of jl± by replacing the 
symbol ji± by cjj,± everywhere on both sides of equations (3.1) and (3.2), for some numerical 
constant c. 

In the preceding equations, Cf and Ca are the quadratic Casimirs, and dp and dA are 
the dimensions, of the fundamental and adjoint color representations. It's also convenient to 
define the trace normalization factor by tr (T]jT|.) = tuS^^, where are color generators. 
In general, tR — dRCu/dA and = Ca- For QCD, 

Cf = I , Ca = 3, c^f = 3, dA — 8, tp = |, ^a = 3. (3.4) 

For QCD with iVf massless fermion flavors, the Debye mass is given by 

ml = {tA + m^) y'T' = (1 + |A^f) g'T'. (3.5) 

Also, the Jg^qq formula in (3.1c) is for a single flavor of quark in the flnal state, and so (1.1) 
should be multiplied by a factor of A'f if one wants the total rate for a gluon to spht into 
any qq pair. 

One can formally expand in powers of inverse logarithms by solving (3.2) by iteration, as 
discussed qualitatively in section Sec. II B. If some initial guess Q±o ~ (P/Ty^'^ is made for 
(1±, then the iterated approximations are 



Ail — -R(Oio)' 

Ai2 = Rifili) = R{R{Qlo)l 



(3.6) 
(3.7) 
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where R{t^±) represents the right-hand side of (3.2). The result pL±_2 is vahd to next-to- 
leading-log order and is actually the form we derive our result in later in this paper. If 
one changes the initial guess Q±_q by an 0(1) multiplicative factor, it only affects the result 
for fi^2 at yet higher order in the inverse logarithm expansion. That's fine as a theoretical 
statement, but it leaves ambiguous how to choose Q±q for numerical comparisons such as 
Figs. 1 and 2. We therefore re-organized our result into the natural form (3.2), which 
provides a specific prescription for determining fi^. The slight cost is that, instead of a 
closed form expression hke (3.7) for the NLL result, we have a simple imphcit algebraic 
equation (3.2) for fij^ that needs to be solved numerically. 



IV. DETAILS OF THE CALCULATION 
A. The equations to solve 

We will work in the formalism of Ref. [9, 10]. In their notation, the functions 7a^6c 
describing nearly coUinear splitting/ joining are 

_ p''^ -\- _ _ ~ 

Iq^qgiP'^P^ k) = lq^qg{p']P, k) = j ^q{p',P, k) , (4.1a) 

- k'^ -\- p^ - 
l9^qqip'\P,k) = ^^^:F^{k,-p,p') , (4.1b) 

p'4 _|_ p4 _|_ ^4 

^g^gg{p';P,k) = -,3-3^3 ^^^ip' , P , k) , (4.1c) 

where 

^s(p',P,k) ^ J -^^2h.ReFs{h;p',p,k). (4.2) 

Here Fs is the solution to the following integral equation 
2h =i 5E{h; p\ p, k) Fs{h; p\ p, k) 

+ 9' J ^>^(«^){^^A mh-p',p,k)-Fs{h + p'qr,p',P,k)] 

+ {Cs - ^Ca) [Fs{h;p',p,k) - Fs{h -kqr,p\P:k)] 

+ ^Ca [Fs{h;p',p,k)-Fs{h-pqr,p',P,k)]^, (4.3) 

which we will refer to as the LPM equation. In this equation, 

SE{h;p.p,k) = — + — - — + — (4.4) 

represents the energy denominator Eg^^ + Eg^p — Es^p' in a ^ pk splitting/joining process. 
Here, rries^s is the 0{gT) effective thermal mass of hard particles of species s.^ The two- 
dimensional vector h is related to transverse momentum and physically corresponds to the 



We will not need explicit formulas for meff,s in this paper. 
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combination 

h = kp± — pk±. (4.5) 

In terms of the qualitative discussion of Sec. II, one can crudely think of the order of 
magnitude of h as representing h ~ PQ±. The function A{q±) is the integrated correlator 



Aiq±) = J '^{{A-iuj,q^,q,)[A-iuj,q^,q,)r)) 



(4.6) 



where A~ = A^ — A^ and ((AA*)) is the thermal Wightman gauge field correlator (neglecting 
the momentum-conserving 5-function) . If momentum transfers are small compared to typical 
plasma particle momenta, then one may evaluate the self-energies in this correlator in the 
hard thermal loop approximation. For the case of equilibrium, A then has the simple form 
[14] 

As discussed in Ref . [9] , the last formula holds more generally for the case of any homogeneous 
plasma where the distribution of plasma particle momenta is isotropic. The "temperature" 
and Debye mass to use in (4.7) in such a situation are^ 

y ^^ S.->AJ^/.(P)11±/.(P)] 

ml = 2g^Es'^stsI^^, (4.9) 

where = i^s/ds is the number of degrees of freedom of a species s excluding color; fsip) is 
the phase space distribution of particles of type s per spin and color degree of freedom; and 
the species sum is over gluons, flavors of quarks, and flavors of anti-quarks. For our analysis 
of splitting processes, this generalization to isotropic non-equilibrium situations will only be 
valid if the plasma particle distribution functions do not significantly change over the course 
of a formation time. 

The simple form (4.7) for A{q±) is justified by our parametric assumption P <^ 
T/fcKg ln(Q;~^)]. We will give some discussion of what happens at higher P in section VI, but 
we leave a full NLL calculation at higher P for future work. 

In (4.3), the variable q± represents the transverse momentum exchange from individual 
2— >2 scattering processes. Solving the integral equation for Fs then accounts for summing 
up multiple scattering into the LPM effect. For a discussion with notation similar to that 
used in this paper, see, for example, Ref. [8]. 



B. The leading log solution 

The Coulomb logarithm arises in the LPM equation (4.3) from the region of q± integration 
where mo <S q± with q± still small enough that h + p'q± , h — kq± and/or h — pq± are 



For QCD with N{ flavors, and identical distributions of all quarks and anti-quarks, this would be = 
9'!^ [6/.(l + /.) + 2Arf /,(1 - f,)]/ml and = 29^J^ [6/, + 2Nd<,]/p. 
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still close to h. Roughly speaking, this is the approximation that the important momentum 
transfers q± from individual collisions are large compared to niD but small compared to the 
total Q± from all the collisions during the formation time. Deep in the LPM regime, there 
will be more and more scatterings in a formation time, and so Q± ^ q±. 

In this same limit, the last term of (4.4) will dominate, and we can approximate 

We can also approximate the differences of l^^'s on the right-hand side of (4.3) by Taylor 
expansions, keeping the first term which does not integrate to zero by parity. The result is 

xWlF,{h;p',p,k) J ^,QlA{q±). (4.11) 

Prom (4.7), the remaining qr^ integral is logarithmically UV divergent. In the original 
integral, this divergence is cut off when our approximation that h + p'q_L etc. are close to h 
breaks down. Let q± ~ Q±o represent any rough estimate of this breakdown scale. Then, 
to leading order in logarithms, 

where Qj_q = Q_lo/^d- The differential equation (4.11) for the leading-log approximation to 
Fs{h) can then be solved, applying the boundary conditions that Fg remain finite at /i = 
and /i — > oo. The result is^° 



Fso(h) = i Ap'pk 



exp ( -e^"-''^^ - 1 



(4.13) 



where 



Hi = \^^mlT\p'pk\ [i^Ap'^ + {Cs - ICa)^^ + \C^p^] In(gio)} ' ■ (4-14) 

The ± sign in (4.13) should be chosen as the sign oip'pk [which is negative in the case of 
(4.1b)], but this sign will not have any effect on our final answers. In terms of our previous 
qualitative discussion in Sec. II and qualitative identification of h ~ PQ±, the above squared 
width of the distribution (4.13) is parametrically of order ~ {PQ±y ~ P^m^Q\^ 
with Ql given by (2.12). 

Plugging (4.13) into Eq. (4.2) for JF,, one finds the leading-log approximation 

J^s^^V2\p'pk\Hl (4.15) 



In solving (4.11), it is convenient to use rotational invariance to first write Fg as h times a scalar function 
of h\ 
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Combining this with the equations (4.1) for the sphtting functions 7 gives the leading-log 
approximation to our result, which just corresponds to replacing all three logarithms on the 
right-hand side of (3.2) by ln(Q5_o)- 



fi'i{xi,X2,Xs; 81,82,83) 
~ mo 



2Hl_ 



"2 


P" 


1/2 




- X1X2X3 








TT 







S3; ^3 



1/2 



in(gio) 



(4.16) 



Our leading-log result is the same as that derived by other authors/^ which were based 
on general formalisms for the LPM effect but applied in a "static" approximation where the 
color fields of plasma particles were treated as simply screened by a Debye mass in the same 
way as static electric fields. In this approximation, the correlator A{q±) of (4.7) is replaced 
by^^ 



«4.static(?±) 



2 \2 ■ 



(4.17) 



The difference between this static approximation and the actual case is that the plasma 
screening of the non-static color electric and magnetic fields generated by moving plasma 
charges is different from the screening of the static, purely electric fields. However, the 
leading log result is generated by q± ^ mo, and in this case all screening effects can be 
ignored. Accordingly, (4.17) is the same as (4.7) when q± ^ rriD. When we go to NLL 
order, individual momentum transfers q± of order rriD will be important, and we will find a 
difference between a full treatment of plasma screening and the static approximation. 



C. The NLL solution 

To discuss the expansion, it is useful to introduce notation 

(/l,/2) = / ^'^^^ ■ ^'^^^ ^^-^^^ 

for the inner product of two vector functions of h. In this notation, the basic LPM equations 
(4.2) and (4.3) become 

and 

S = CF, (4.20) 



In particular, our leading-log result for dVg^gg/dx using (1.1), (3.1a) and (4.16) is the same as Eqs. 
(19-20) of Ref. [11] with their \n{{k\)/ml,) identified as our hi{Q\^). 

See, for example, Eq. (2.17) of Ref. [3]. The difference by an overall factor of ttT is because they have 
normalized their version V{q\) of A so that / (Pqi_ V{q\) = 1. See also Ref. [15]. 
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where we've defined the "source" S of the last equation by 

S = 2h (4.21) 

and C is a linear operator defined such that Ci*" is the right-hand side of the LPM equation 
(4.3). 

Let Co represent the approximation to C which gives the leading-log result. Expanding 
to first order in powers of 5C = C — Co, we can write 

(S,F) = (S,C-iS) 

= 2(S,Co-iS)-(S,Co-^CCo-^S) 

= 2(S,Fo)-(Fo,CFo), (4.22) 

and so 

Re(S, F) = 2Re(S, Fq) - Re(Fo , CFq). (4.23) 
The first term is just twice the leading-log result, which corresponds to 



Re{S,Fo) = ^\p'pk\H' 

,2 (I 



TT 



^ \p'pkf' { [|Cap'^ + (C, - \C^)k' + |CaP^] In(gio)}'^' 



(4.24) 



We now need to compute 

Re(Fo , CFo) = Re/i + |CARe72(p") + (C, - |CA)Re72(A;') + ICpJiehip") , (4.25) 

where ^ 

Re/i = j ^ Re[Fo(/i) ■ 1 6E{h) Fo{h)] (4.26) 

and 

^^^'^'^ j W?^' "^^^^^ ^"^^^ ■ ^^"^^^ ~ ^"^^ ^ ^"^"-^^ ■ ^^'^^^ 

In doing the Ji integral, one can ignore the tv? terms in Eq. (4.4) for 5E and use the 
approximation (4.10). The w? terms in (4.4) are suppressed by order w?P'^/h^ compared 
to the /i^ term. In terms of the earlier qualitative discussion of section II, h ~ PQ±, and 
so, from (2.13), this suppression factor is of order 



m-^P^ 1 



T \T 



-1/2 



(4.28) 



It is suppressed by a power of P/T, not simply a power of the logarithm ln(P/T), and so 
these mass terms do not contribute at any finite order in our inverse log expansion. The Ii 
integral is trivial and then gives 

Re/: = ^^=|Re(S.f.). (4.29) 

TTV 2 ^ 
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The evaluation of I2 is briefly outlined in the appendix and gives 
where 



- e^-/^ ^ . (4.31) 

Parametrically, is the same order as (4.28) and so is small compared to one. We therefore 
only need the small u,^ expansion of the integral (4.30)/^ which is derived in the appendix: 

= (2 - 7e - ln«.)«. + 0{ul), (4.32) 

TT 

and so 

= -^^^^ (2 - 7e + ^ - In W\) W\ + 0{ul). (4.33) 
Using (4.14) and (4.24), this can be rewritten as 

Re72(«') = -- Re(S, Fq) ^ ^"^^^ ^ . , (4.34) 

where ^ is the NLL constant under the log defined in (3.3). Putting our results for Ii and 
I2 into (4.23) and (4.25), the NLL result for Re(S,F) is 

Re(5, F,) = Re(5, Fq) 

2 I iCAP'^ ln(Qi„) + (C, - iCA)*:2 ln(Q2 ^) + ICaP^ MQlo) 

with i^s given by (4.14). At leading-log order, where all logarithms are treated as the same, 
the multiplicative correction above simply reduces to a factor of one. At NLL order, its 
effect on the leading-log result (4.24) is simply to replace the curly brackets in (4.24) by 
the numerator of the big fraction in (4.35). This corresponds to the iterated NLL solution 
(3.7) quoted in section 111, when the splitting functions 7 are written in the form of (3.1). 
Note that in (4.35) depends on Q±o- As mentioned earlier, the self-consistent equation 
(3.2) provides an answer which is equivalent at NLL order but which does not depend on 
specifying an initial guess Q±_q. 



D. Dynamic vs. static treatment of screening 

Since the static approximation (4.17) to screening is sometimes used in the literature, we 
will take a moment to discuss how the NLL order result would have been different had we 



The exact result for (4.30) can be expressed (somewhat uselessly) in terms of the Meijer G function as 
h{n^) = k + i (7e + \nu.? - (n« \,'^],) 
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FIG. 4: Similar to Fig. 2 but for ^q^gq. x is the momentum fraction of the gluon. 
made that approximation. This is simple to do by comparing with (4.7) and noting that 

d{mf)} d\u^\ 

where the last equality uses (4.31) considering H as fixed. The only change comes in the 
calculation of J2. Applying (4.36) to (4.33), 

RehiK^) - Relt'-^iK^) = -^^^^ (l - 7e + ^ - In K\) K\ + 0{ul). (4.37) 

This produces the same NLL result (3.2) as the fully dynamic case but with 

e ^ U.tic = exp (1 - 7e + f ) - 3.3474 . (4.38) 



V. NUMERICAL RESULTS 

In the introduction, we have already discussed for g gg the comparison of NLL results 
with a full computation at leading order in powers of ag. The scaling in Fig. 2 was chosen by 
inspection of the small x limit of the NLL result given by (3.1) and (3.2), choosing powers 
of X so that the curves will be similar for small enough x. We give similar results for "jq-^qg 
and 'Jg^qq in Figs. 4 and 5. The fact that the small x curves in these figures do not fall 
closer together is because we have not quite gone to small enough x and because we have 
not bothered with factors of ln(l/a;) in our consideration of how to scale the axis. 

For the case of g — gq, there is no final state symmetry that relates results for x to 1—x, 
so we show X > 0.5 results for this case in Fig. 6, scaled so that the curves are similar in the 
1—x — i> limit. 

From these graphs, one can confirm for all splitting processes the general claim made in 
the introduction that the deviation of the NLL approximation from the full result at leading 
order in powers of as is roughly 20% or better for p< > lOT. 
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FIG. 5: Similar to Fig. 2 but for ^g^qq- 
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FIG. 6: The case q ^ gq iov x > 0.5. 



VI. THE CASE OF LARGE ln(P/T) 



Throughout, we have assumed that as is small enough that we can ignore all effects 
suppressed by powers of as. In particular, we have ignored the issue of what renormalization 
scale as should be evaluated at: m^, T, Q^, or P? Formally, 



asin') = asi/J,) 



/2 



l+/5oasln( ^ ) +Oiai\n') 



(6.1) 



17 



and so we can ignore the ambiguity of scale choice if ccg ln(//'/ /i) is small for relevant different 
possibilities fi and fi'. 

In this paper, we have studied the deep LPM regime, where lia{Q±/mB) ~ ln(P/T) is 
large. Nonetheless, we have so far implicitly assumed that ttg ln((5±/mD) ~ asln{P/T) is 
small, since we have ignored all higher order corrections in as, including the running of 
the coupling constant. This assumption is also implicit in the full leading-order results of 
previous work, such as the LPM equation (4.3). However, in applications there could be a 
significant difference between as{m-D) and as{Q±) — that is, asln{Q±/mY)) may be large. 
In this section, we discuss what to do in that situation. 

Earlier, we made the parametric assumption P -C T/[Q;g ln(Q;~^)] in order to simplify the 
analysis of this paper, which is a first step towards computing the NLL result for the general 
case. With this assumption, the expansion parameter asln{P/T) is never large in the weak 
couphng hmit, since then asln.{P/T) < asln{a~'^) <S 1. In general, different choices of 
scale fi in as(//) are formally not significant in the weak coupling limit if the different scales 
only differ by some power of the coupling. So, for example, there is no significant difference 
between as{m-£,) and as{T) in the weak coupling limit. 

In this section, we will abandon the restriction that P <S T/[agln{a~^)] and consider 
extremely large P for which as{Q±) may be significantly different than Q;s(mD) — Ois{T). 
We will correspondingly restrict ourself to a leading-log analysis, though now resumming 
the large factors of ctg ln((5±/mD) due to running of the coupling constant. 



A. Large q± behavior of ^(gj.) 

In our previous analysis, the leading- log result was dominated by uib <S q±, where the 
gauge field correlation function had the simple form 

Trn^ 

A{q±) ~ (mo < ?± < T). (6.2) 

?± 

This can be understood as arising from the lower part of the scattering diagram of Fig. 
3. The l/q^ is the contribution (— a;^ -|- q^)~^ of the gauge propagator to the rate, with uj 
set equal to as in the general definition (4.6) of A{q±). Setting a; = is a reflection of 
energy conservation for the high-energy particle (the top solid line of Fig. 3) in the limit that 
Q <^ P. The factor of Tm^ in (6.2) reproduces the contribution of the plasma particle and 
its phase space integral, integrated over q^ as in (4.6). It's appearance is more transparent 
if we write an explicit formula for the Debye mass, 

, . r ^3/) 

Tml ^g'J2 ^stR. J ^ fs{i) [1 ± fs{i)]. (6.3) 

Here one can see the g'^ associated with the contribution to the rate of the bottom vertex 
in Fig. 3, the Bose or Fermi distribution f{£) for the probability of encountering a plasma 
particle with momentum i, a final state Bose enhancement or Fermi blocking factor 1 ± 
f{\£- + q\) — l±/(£), a sum over species s of the plasma particle, and the appropriate group 
factors and numbers of degrees of freedom. 

For the opposite limit of q± ^ T, one obtains a similar expression, but there is no final 
state factor 1 ± /. That's because if £ ~ T is a plasma particle momentum, and q± ^ T, 
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then |£ + qj > T and so 1 ± f{£ + q) ~ 1. The form of A{q±) in this hmit is 

A{q^)^^ (6.4) 

where 

^^Y.-stR.j J^Jsii) (6.5) 

is a measurement of the density of plasma particles, weighted by group factors. For QCD 
with A^f massless fermion flavors, N is 

= if ('^ + ^"^^'0 = H + -'""'^ ^^-^^ 

where C,{s) is the Riemann zeta function. Compare this formula to (3.5) to see the difference 
between Tm^ in (6.2) and g^M in (6.4). 

For 3-flavor QCD, g'^M is about 15% smaller than mpT". For any phenomenologically 
relevant values of the coupling ccg, this difference is unlikely to be significant compared to 
other corrections, such as higher-order effects. So, one could reasonably just start from our 
earlier results based on the original formula (4.7) for A{T). However, the conceit of this 
paper is to work out precise results in the formal limit of arbitrarily weak coupling (and 
large logarithms), and so for sufficiently high energy jets we use (6.4). 

Now consider our previous leading-log analysis. The Coulomb logarithm was generated 
from (4.12): 



/ 



^"^ A AM n ^ , M „ I ^ I (e, « D. (6.7) 



For Q± ^ T, there are instead two logarithmic contributions, coming from the different 
integration regions represented by (6.2) and (6.4): 



27r V"^D / 27r \m-D 

If P ^ T/[Q;g ln(Q;~^)], so that Q±/T 3> T/niD by (2.9), then the second logarithm domi- 
nates. Since in this section wc will be interested in P large enough that there is significant 
running of the coupling constant (and so formally P larger than T times any fixed power of 
l/cKg), we shall henceforth assume this is the case. If we ignored the running of the coupling, 
we'd have 



/ 



(271)2-^^ v-^-^ J^^^ 27r ql 27r V^d^ »tHo^s') 



in this case. Given our assumpations and approximations, we could have just as well made 
the lower limit T instead of mo in (6.9). However, because of the small difference between 
g'^Af and m^T, taking the lower hmit to be mo is probably shghtly better in practice. 
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B. Running of coupling with q± 



We begin with a discussion of the scale of the coupling g in the g^A{q±;g^) factor in 
the LPM equation (4.3), where we now write A{q±:, g'^) instead of A{q±) to emphasize the 
fact that A depends on g'^, e.g. as in (6.4). Physically, this factor is proportional to the 
rate for a high-energy particle to scatter off of a plasma particle, with momentum transfer 
q±. (See, for example, the discussion in Ref. [8].) A momentum transfer of q± corresponds 
to an impact parameter between the high-energy particle and the plasma particle of order 
l/q±. If this distance is very, very small, then g'^ in this Coulomb scattering amplitude 
should be correspondingly smaller because of the anti-screening of the QCD vacuum. (For 
^ the Debye screening length 1/mD, medium effects are ignorable, and so the relevant 
consideration for such a collision is of the running coupling constant in vacuum.) The 
upshot is that the appropriate scale for g"^ in this scattering rate should be set by q±. This 
is precisely what we would get by summing up all 1-loop bubbles on the exchanged gluon 
line. We therefore expect that the LPM equation (4.3) should be modified by replacing g^ 
by g^{q±) so that it becomes 

/j2 
j^9'{q±)A{q±;g'{q^))[---]. (6.10) 

This is the same prescription as used in a related context by Peshier in an analysis of 
coUisional energy loss [16]. 

Now repeat the leading-log analysis we reviewed in Section IV B. Correspondingly, we'll 
use the 1-loop renormalization group result for g'^{q±), 



-(3o ln(//VA2) 



(6.11) 



For QCD, A represents Aqcd, and 



- (IICa - 4iVftF) (33-2iVf) ^ 

Po — rrr-^ — TTT-^ < U. 



487r2 



487r2 



(6.12) 



(For QED, A represents the ultraviolet Landau pole, and Pq > 0.) The integral analogous 
to g^ times (6.9) is 



/ 



(Pq±_ 



g\q±)qlA{q^;g\q±)) ~ 



q^ dq^ g\q^)N 

ql 

dq± 





27r 




r^Q± 




~m,D 


M \ 





87r/?2 



q^\n\qJK) 
1 1 



,ln(mD/A) ln(g^/A), 



(6.13) 



Following through the previous leading-log derivation, we will get the same result for Fq and 
H but with the replacement 



g'miT\n{Qi,) 



ln(mD/A) ln(gxo/A), 



^ [9\m^) -y (Q±o)] ^g_-L4^ 
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in (4.14). The correspondence with the original leading-log result is easier to see if we use 
(6.11) to rewrite this in the equivalent form 

/m^rin(gio) ^ g\m^)g\QA)MHQl^). (6.15) 

It is interesting to note that, because of asymptotic freedom in QCD, the running-coupling 
formula on the right-hand side of (6.14) is finite as P (and so Qj,) becomes infinite. In this 
limit, the value of ois{Q±) is irrelevant — a fact slightly obscured by (6.15) but made clear 
by (6.14). In this case, the momentum transfers in individual 2^2 scatterings which 
dominate the integral (6.13) range, roughly speaking, from T up to those where as{q±) can 
first be considered small compared to as{T). In particular, once Q±_ is large enough, the 
scale of g_L does not continue to grow as one increases P and (5±. 



C. Remaining g"^ and synthesis 



There remains one other factor of in the problem, which is the cost for emitting the 
energetic, bremsstrahlung particle. Because of the LPM effect, a bremsstrahlung gluon 
cannot resolve individual 2— >2 collisions with the plasma but is only sensitive to the net 
deflection of the high energy particles over the entire formation time. The relevant scale for 
the g'^ cost of bremsstrahlung is therefore plausibly the total momentum transfer^"^ (2-13) 
{P/Tf/Hn^'\P/T). [We will ignore the issue of whether this last P should be P 
or xP or {1 — x)P but will assume that x and 1 — x are large enough that there is not much 
difference in ccg.] We therefore propose that the correct leading- log formula in the case of 
small Qfs but large «sln(P/T) is given by (i) replacing ctg by as{Q±) — Q;s(mDQ±o) in (3.1), 
and then (ii) using (6.14) to replace (4.16) by 



1/2 



2 P 

- XiX2X^ — 
TT 1 



1/2 



X 



X 



1/2 



1/2 



(6.16) 



[Note that it is the combination /i^ = m^fi^ which appears in (3.1).] Here, Q±o is chosen 
with order of magnitude given by (2.13), or one could simply self-consistently choose Q^o — 

For the case of large g"^ ln((55_o)) there is not an obvious benefit to pushing further to find 
the NLL solution to the running LPM equation (6.10). For large g^ln^Q^g), the expansion 



1/Q± is the transverse distance corresponding to the deflection of the high-energy particle due to 2^2 
scatterings during a formation time. To see this, consider that if the particle picks up transverse momen- 
tum Q± in that time, then it's transverse velocity will be of order Q±/P. Multiplying this transverse 

velocity by a formation time fform ^ P/Q\ (2-4) then gives 1/Q±. The same scale is also the scale of 
the quantum mechanical uncertainty of the transverse position during a formation time (which is why 
LPM interference can occur). In contrast, the scale l/q± relevant to section VI B is the typical transverse 
distance between the high energy particle and the plasma particles it is scattering from. 
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parameter \\.'n{Q±/m-D)]~^ is parametrically the same order as as, and so one cannot obvi- 
ously justify neglecting corrections that are higher order in as that have not been included 
in the running of the coupling. On the other hand, it would be nice to have a single, well- 
defined formula that interpolated our previous NLL result into the realm of large g'^ ln((55_o)- 
But we shall not pursue this here. 

We have made plausibility arguments about which renormalization scales should be used 
for evaluating the coupling g^. It's possible, however, that if one starts looking at corrections 
suppressed by ag that there are other large logarithms which arise, unrelated to the runnings 
we have described. To settle the issue definitely, it would be nice to have explicit weak- 
coupling calculations beyond leading order in ctg of a well-defined physical quantity that is 
dominated by particle splitting processes. 
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APPENDIX A: I2 
1. The I2 integral 

Here is one way to do the integral hii^"^) of (4.27). First use the standard trick of rewriting 

-= / dAe-^(«'+-'), (Al) 

so that 



_ g-A(<zi+m2,) 



Next rewrite (4.13) in the form 



Fo{h)-[Fo{h)-Fo{h + Kq^)]. 

(A2) 



(A3) 



where C = iAp'pk, A = e^^'^^^/H'^, and the limit e — > O"*" is taken at the end of the day. By 
introducing e, we can split I2 into 



(A4) 
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where 

I{^)=9'TjyX{l-e-'<) J-^2j0^2'-"''^^'Fo{h).Fo{h + j) (A5) 

and we have introduced the notation j = —Kq±. The integral I would be divergent if we 
had not introduced some sort of regulator like the e. Now note that (A5) has the form of a 
convolution, and so by Fourier transformation we can recast it as a single two-dimensional 
integral over a Fourier conjugate variable which we'll call B. Two-dimensional Fourier 
transformation of h or j takes 



^2 



and so the integral becomes 
By (A4), then 

= ^(1 - / ^ (1 - e-'^^^^/^^) fe-^^/^^ - e'^^/^^V . (A9) 

We can now safely take the limit e ^ 0, giving 

Finally, doing the B integral and changing variables from A to r = 2\/Ak'^ gives (4.30). 
2. Small Uk, expansion 

To expand (4.30) in powers of it is convenient to notice that /2(0) = and then 
instead expand 

(All) 







where Ei is the exponential integral. The small expansion is 

^(l-7E-ln««). (A12) 

Integrating both sides and using /2(0) = then yields (4.32). 
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